Magnetic anisotropy of BaCu2Si207: theory and antiferromagnetic resonance 
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Antiferromagnetic resonance (AFMR) of BaCu2Si207 and a microscopic theory for the magnetic 
anisotropy of spin 1/2 chain compounds with folded CuOs geometry being in good agreement with 
the available data are presented. The AFMR studies at 4.2 K show the existence of two gaps (40 and 
76 GHz) at zero magnetic field and of two spin re-orientation transitions for H || c. The microscopic 
origin of the two gaps is shown to be Hund's rule coupling which leads to a "residual anisotropy" 
beyond the compensation of the Dzyaloshinskii-Moriya term by the symmetric anisotropy which 
would be valid without Hund's coupling. 
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There exist several cuprate compounds with CuOa 
corner-sharing chains that are known as nearly ideal 
model compounds for onc-dimensional spin 1/2 systems. 
Those with folded CuOs zigzag chains, like BaCu2Ge207 
(Neel temperature of Tat = 8.8 K) |l| or BaCu2Si207 
(Tat = 9.2 K) |, |, |, I open the fascinating possibil- 
ity to study noncoUinear magnetism or the influence of 
a Dzyaloshinskii-Moriya (DM) term in a quantum spin 
chain. Especially BaCu2Si2 07 attracted much interest 
recently in connection with the finding of two consecu- 
tive spin re-orientation transitions for a magnetic field 
applied along the easy c-axis, which were discussed con- 
troversially, however jj, ||. Furthermore, one or two 
gaps 1^ in the spin-wave spectrum measured by neutron 
scattering were reported. And also the microscopic ori- 
gin of the magnetic interaction energies was not clarified 
up to now. The present work addresses all these points 
by a combined experimental and theoretical study, us- 
ing antiferromagnetic resonance (AFMR) measurements 
and perturbation theory with respect to spin-orbit (SO) 
coupling in a Cu-O-Cu cluster with folded geometry. 

Our experimental study shows clearly the existence of 
two different gaps at zero magnetic field in BaCu2Si207 
and a magnetic structure evolution with increasing field 
H \\ c with two re-orientation transitions corresponding 
to a spin rotation towards the middle 6-axis followed by 
a rotation towards the hard a-axis which confirms the 
neutron scattering structure of Ref . Q . 

Theoretically we show that the deviation from a 
straight bonds in a Cu-O-Cu cluster leads to a DM inter- 
action and a spin canting of neighboring spins at T < Tjv . 
The canting angle is rather large and the DM interaction 
alone would lead to much larger gaps in spin wave spec- 
trum and AFMR than actually observed. It is the sym- 
metric anisotropy part which compensates the DM one, 
as proposed in Refs. ||^, |7| and applied here to the given 
situation. A clear experimental confirmation of this effect 
was found in Ba2CuGe207 H. We argue that Hund's 



coupling leads to a "residual anisotropy", beyond this 
compensation, which agrees well with the experimental 
data and explains the observed easy axis behavior and the 
two gaps mentioned above. Naturally, our theory is also 
valid for BaCu2Ge207 which differs from BaCu2Si207 by 
the interchain couplings. The more simple AFM struc- 
ture of the Ge-compound allows for a weak ferromagnetic 
moment which gives direct information on the spin cant- 
ing angle 0]. 

The crystal structure of BaCu2Si207 and 
BaCu2Ge207 belonging to the orthorhombic space 
group Pnma is made of almost isolated CuOa corner- 
sharing chains running along the c-axis Q . The Cu-O-Cu 
bond angle is found to be 124° in the case of BaCu2Si207 
which is smaller than that of BaCu2Ge207 (135°). For 
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FIG. 1: The frequency vs field diagram [H \\ c) of AFMR 
modes in BaCu2Si207 at T = 4.2 K. Dashed lines: AFMR 
modes of a biaxial antiferromagnet without the DM interac- 
tion in a coUinear (phase 1) and a canted (phase 2) phase 
according to Ref. [[ll[ . 
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both compounds a broad maximum was observed in the 
temperature dependence of the magnetic susceptibihty 
indicating ID magnetic behavior |^ and allowing to 
extract the interchain exchange constants J — A7 meV 
(J = 24 meV) for the Ge (Si) compound. The single 
crystals of BaCu2Si207 used in an AFMR experiment 
were grown by a floating-zone method. Since the a- and 
c- axis lengths are almost the same, particular attention 
was paid to the crystal orientation. The AFMR study of 
BaCu2Si207 was done using a simple millimeter-range 
video spectrometer The microwave sources were 

either Gunn's diodes or back-ward wave tubes. The 
field dependence of AFMR modes in BaCu2Si2 07 at 
T = 4.2 K is shown in Fig. 1. Two gaps are clearly 
seen at H — 0, vlf = 40 GHz (low-frequency mode) 
and vhf = 76 GHz (high-frequency mode). The gaps 
values are about 10% smaller as compared to the earlier 
reported values 0.21 meV and 0.36 meV measured at 
lower temperature T = 2 K ||^. The observed frequency 
softening of the low-frequency mode at Hd and Hc2 
confirms the existence, for this particular magnetic field 
orientation, of two consecutive spin re-orientation transi- 
tions. As it was first proposed in Ref. this behavior 
is characteristic of a AFM with the DM interaction in a 
magnetic field applied along the easy axis and the DM 
vector orientation along the middle axis (c- and &-axis 
respectively in the case of BaCu2Si207; see Fig. 4 of 
Ref. pO| for comparison). Furthermore, our observation 
of quadratic field dependence of the high-frequency 
AFMR mode for H \\ a (not presented in Fig. 1) allows 
the conclusion that the a-axis is the hard axis of the 
magnetic anisotropy tensor of BaCu2Si207 p^ . 

To understand the microscopic origin of the observed 
behavior in Fig. 1 and also other experimental facts al- 
ready established for the two compounds in question we 
investigate the special Cu-O-Cu bond sketched in Fig. 
2 with an arbitrary bond angle it — (j). The isotropic 
exchange, the DM term and the symmetric anisotropy 
are calculated by perturbation theory with respect to 
SO coupling and kinetic energy. The Hamiltonian con- 
tains all the 3d orbitals at the 2 Cu-sites (denoted by m 

or m' = {0{d^2_y2),z{d3:y),x{d3:z),yid^z)Aid3z^-r^)}, 
where the usual notation is given in parenthesis) and the 
three 2p orbitals {n or n' = {px-,Py,Pz}) at the oxygen 
site in between. The local, unperturbed part contains 
the Coulomb terms 



Ho = EPnP + ^{nPnP 
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with 



~ X/m ^ctm ~ X/mo- d^imcrdama and vP — 

and where a = A ov B (sort of Cu), and 
cr (cr') are spin indices. The charge transfer energy Ap = 
EP - Eq has been set to 4 eV |3) and E^ - Eq = Ed = 2 
eV shall be the same for all m = {x,y,z,l} |H. We 




FIG. 2: The orbitals of the CU2O cluster and the coordinate 
systems used in the calculation. The two ground state orbitals 
Sd^2_y2 and 3d^,2_yr2 are fixed by the positions of the two 
CUO4 plaquettes perpendicular to the plane of the Figure. 



use in the following the standard parameters of cuprates, 
i.e. Up = 4 eV and Ud = 10 eV. All the remaining terms 
of the Hamiltonian are treated as a perturbation. That 
concerns especially the kinetic energy 



t = ^ tam,n{d^amaPn<y + h.C.) 



(2) 



and the spin-orbit interaction 



A 



Hso ~ 2 ^ ^ ^ . ^am' cr'-^mm'^CTCT' 'Earner , (3) 
ha mm'aa' 

(A = 0.1 eV) with the only nonzero matrix elements 
L^o = 2z, = -z, and ^ -z {L^^, = {L^.J*), 
where 6 is a Cartesian coordinate and cr^g./ the corre- 
sponding Pauli spin matrix. The Hund's exchange inter- 
action is given by 



aa amm' 

-JhpY^ E PiaPna'pi,^,Pn'a 
era' nn' 



(4) 



In a first step, the intermediate oxygen orbitals are ex- 
cluded and an effective spin Hamiltonian 



2 E ^!lmicr'^'4m3(T'rfB„i2CT'^-B™4CT 
{m}aa' 

■J{miA,m2B;m3A, ruiB) 



(5) 



is derived. The necessary exchange terms are calculated 
in 4th order perturbation theory with respect to the ki- 
netic energy Ht. The magnetic anisotropy terms arise 
due to the SO interaction. Up to second order, the result 
can be written as 

Hdm + Ha = D{Sa X Sb) + Sa^Sb , (6) 

with Vl"'' = F'^'' + F'"' - (^"''(Ec ^'"')- The first order term 
corresponds to a DM interaction: 

zAL|5„j(J(0A,mB;0A,0B) - J(mA, Ob; 0^, Os)) 



E„i — En 



(7) 
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and the second order describes a symmetric spin 
anisotropy 
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LpmL'^n'oJ Ob; m'^,OB) 
{Em — Eo){E„i' — Eo) 



EornLpm' J{mA,m'g;OA,OB) 

{Em ~ Eo){Em' — Eq) 
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There we neglected those terms in second order ^ 
which lead only to a correction of the isotropic exchange 
interaction. To apply the general theory to the given sit- 
uation one has to specify the transfer terms that couple 
\A, 0) and \B, 0) due to the different perturbations. First, 
we neglect Hund's exchange. The oxygen p^, Py and pz 
orbitals are defined in the local coordinate system of Cvla 
{x-y-z in difference to the x'-y'-z' system at Cus) and 
the relevant transfer terms between Gwa and the inter- 
mediate oxygen are tAo,p^ = tpd and tAy,p, = tpd/V^, 
where we used the ratio between a and tt transfer to be 
1 : \/3 which is valid with good approximation. To obtain 
the transfer integrals of oxygen with the Cus orbitals we 
have to rotate the \p'^) and \py) orbitals to \px) and \py). 
Then we get the transfer terms 



iso.p^ = cos (f) tpd 
tsy.p^ = sin0 ^ 



tso.p^ = — sin0 tpd 

tpd 

71' 



tBy,p^ = cos ( 



(8) 



The next step consists in calculating the exchange inte- 
grals in H'^^. Let us start with the isotropic exchange 
J = J{Oa, 0_b; 0^, Ob). We use perturbation theory with 
respect to Ht in 4th order. There are two possible in- 
termediate states, the doubly occupied oxygen state \px) 
with an energy 2Ap -I- Up and the doubly occupied cop- 
per state with the energy Ud- Collecting all the possible 
paths we get 



J = AGb"^ with G = 



1 

p 



1 

U'd 



2A„ + Ur, 



(9) 



and b — t^^ cos 0. With the standard value for cuprates 
tpd — 1.3 eV we obtain for BaCu2Ge207 with (/> = 45° the 
estimate J « 95 meV being roughly two times larger than 
the experimental value. The reason is most probably the 
insufficiency of 4th order perturbation theory which is 
also known for the standard CUO2 plane. For the Si- 
compound {4> = 56°) we estimate 60 meV, whereas the 
experimental value is 24 meV. Now, the DM vector D 
shall be calculated. The two ground states |A,0) and 
\B,Q) are coupled by SO and Ht only via the \B,y) or 
the |A, y) orbitals. Therefore, only is different from 
zero and the corresponding exchange paths (in 4th order) 
are shown in Fig. 3 with 



-'^(0a,2/b;0a,0b) 



4Gtp^ cos sin < 



(10) 
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FIG. 3; Exchange paths for the DM term, contributions to 
J(i/a,0s;0a,0s) (a) and J(Oa, Oa, Ob) (b). The initial 
state is indicated by full arrows, the final state by dashed 
arrows. 



being equal to — J(yA, Ob; Oa, Ob). We find 

= 8G6c where c = — sin (/>-^ 
£d V3 

The corresponding canting angle 
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(11) 



(12) 



results in 0.058 for BaCu2Ge207 whereas 0.033 was esti- 
mated on the basis of the magnetization data [Q . 

The contribution in second order of the SO coupling 
(^) corresponds to a symmetric anisotropy term. The two 
relevant nonzero exchange terms are J{yA, Ob; Ua, Ob) = 
— J(yA, ys; Oa, Ob)- The corresponding exchange paths 
are very similar to those shown in Fig. 3. So, the first 
term is characterized by an exchange between \A, y) and 
|i?,0) via the \pz) orbital. The second, mixed exchange 
goes from \A, y) via \pz) to 0) and from \B, y) via \px) 
to I A, 0). The result can be given as 



J(2/a,Ob;2/a,Ob) = 4G 



^/3 



Summing both contributions we get 



(13) 



(14) 



Now, we see that we can combine the isotropic exchange 
(^) , the DM term (^ij) and the symmetric anisotropy ( p^ ) 

4G ((fe2 - ^)SaSb + 2bc{SA X SB)y + 20^5^5^^ 

= JSaSb (15) 

to an isotropic exchange of canted spins S a/b with the 
canting angle 8 (111) Esl. In the zigzag chain one 
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can subsequently rotate all the spins. As a consequence, 
we would obtain a Hamiltonian which is rotationally in- 
variant in spin space and no crystallographic direction 
would be preferred. That is exactly the compensation 
as proposed in Refs. ||, Q]. Now, we are going to show 
that Hund's exchange interaction (|^) leads to a "resid- 
ual anisotropy" . These symmetric anisotropy terms ex- 
ist already for = 0. The anisotropy terms ( pi] ) and 
( p^ ) can be eliminated by the rotation (|l^) and only the 
additional contributions due to are considered fur- 
ther on. Starting with Hund's exchange at oxygen J Hp, 
new exchange integrals become possible in 4th order of 
Ht and first order of Hh- For example, the exchange 
paths for J^'(z^, 0^; z^, Os) are realized due to the trans- 
fer term tj^^.p (= tAy.p^)- It turns out that it is identical 
to JP{yA,OB;yA,QB)^ J'^^: 



J(p) 



J Hp 



A2 

p 



(2Ap 



Up)' 



(16) 



We have no contribution via or \B,x), and there- 

fore, = 0. Next, we consider Hund's exchange at 
copper JhcI- Then the transfer terms tBy,p^ = tBi,p^ = 
cos(j)tpd/V^ become important. One obtains 



j(d)_ 4 ^^rfCos- 
9 A2 



J, 



Hd 



Both terms ( |16| , |17| ) together result in the 
anisotropy" contribution 



(17) 



"residual 
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A2 



J(p) 



(18) 



of one plaquette. To connect the cluster of Fig. 2 
with the real crystallographic structure (see Fig. 1 
in Q), one has to rotate the coordinate axes Cz — 
{za,Zb,Zc) = (0.84,0.38,0.40) and Cy = {ya,-yb,-yc) = 
(0.48, —0.85, —0.21) (and also e^/ and Cy') to the crystal- 
lographic ones. We obtain the diagonal elements 



Yvy 



vl) 



(19) 



with V = {a, 6, c} and only one non-zero off-diagonal ele- 
ment F'^''. With the characteristic parameters Jhp = 0.4 
eV and Jud = 1 eV and setting F^^ to zero, we get 
F'*" = -2 /zeV and F'''' = -0.5 ^eV. These numbers are 
extremely small, but we have checked that they are still 
larger than the classical dipole-dipole interaction. For an 
antifcrromagnetic arrangements of spins along the chain, 
the preferred spin direction is the c-direction which is in- 
deed the case for BaCuaGeaOy § and BaCuaSiaOr [|. 
From the theoretical analysis follows that we can expect 
for BaCu2Si207 an easy c-axis, a middle &-axis and a hard 
a-axis. The diagonal anisotropy energies ( p^ are directly 
related to the two gaps seen in Fig. 1 (and also in neu- 
tron scattering 0]) according to the standard formulas 



vlp=^/2HaiHe and vhf=\/'^Ha2He with He = zJS 
(S - spin, z - number of nearest neighbors) and Hai = 
2z|r''''|5', Ha2 = 22|F""|S' |6). This yields 0.44 meV and 
0.22 meV in reasonable agreement with the experimental 
values 0.36 meV and 0.21 meV. 

In conclusion, our theory is applicable to folded CuOa 
chains and gives the isotropic exchange, the DM term 
and the spin canting angle, the compensating symmetric 
exchange and the "residual anisotropy" due to Hund's 
coupling in good accord with the available experimen- 
tal data of BaCu2Ge207 and BaCu2Si207 and especially 
with the two gaps measured by AFMR for the latter com- 
pound. Possible differences between the two compounds 
might be caused by different interchain couplings which 
were not included in the present study. 
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(PST.CLG.976416) and DFG (UKR 113/49/0-1) for sup- 
port. 
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